Voltage and dephasing probes: a full counting statistics discussion 
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Voltage and dephasing probes introduce incoherent inelastic and incoherent quasi-elastic scatter- 
ing into a coherent mesoscopic conductor. We discuss in detail the concepts of voltage and dephasing 
probes and develop a full counting statistics approach to investigate their effect on the transport 
statistics. The formalism is applied to several experimentally relevant examples. A comparison of 
different probe models and with procedures like phase averaging over an appropriate phase distribu- 
tion shows that there is a perfect equivalence between the models for the case of one single-channel 
probe. Interestingly, the appropriate phase distribution function is found to be uniform. A uniform 
distribution is provided by a chaotic cavity with a long dwell time. The dwell time of a chaotic 
cavity plays a role similar to the charge response time of a voltage or dephasing probe. For multi- 
channel or multiple probes the transport statistics of voltage and dephasing probes differs and the 
equivalence with phase averaging is similarly lost. 



I. INTRODUCTION 

In this work we are interested in the effect of voltagei 
or dephasing probes^ on the transport statistics^ of meso- 
scopic few channel conductors. Voltage probes coupled 
to a coherent conductor permit incoherent inelastic pro- 
cesses: an electron entering the probe is thermalized 
by dissipation and later on fed back into the system. 
In contrast dephasing probes are quasi- elastic and non- 
dissipative, each electron exiting into the probe is reemit- 
ted within the same energy interval. Voltage probes can 
be used to describe a transition^ from purely coherent 
transport to classical inelastic transport whereas dephas- 
ing probes can be used to describe a transition^ from 
coherent transport to classical quasi-elastic transport. 

The main motivation of this work is to provide an im- 
proved understanding of the physics of probe models. We 
achieve this by discussing the models on the level of full 
counting statistics (FCS) and by comparing them with 
each other and in addition with other procedures like 
phase averagingi^ The full counting statistics^ is charac- 
terized by a generating function from which conductance 
and noise and all current cumulants can be obtained sim- 
ply by taking derivatives. For brief reviews on FCS see 
articles in Ref . 0. 

Voltage and dephasing probe models are widely used, 
most often to address the disappearance of quantum ef- 
fects in conductance. The transition from quantum com- 
bination of scatterer to the series addition of resistances 
is the subject of Ref. Q. Other applications include the 
role of equilibration of edge channel populations for the 
quantization of the Hall conductance^ and the transi- 
tion from the quantum Hall effect to the classical Hall 
effect'*, the role of dephasing in the quantum measure- 
ment process^ or the effect of inelastic processes in quan- 
tum pumpingifliii. The success of this approach is due 
to the fact that it can be directly incorporated in the for- 
malism used to discuss the purely coherent limit. This is 
particularly evident in random matrix theor y ^^d^ Since 



the scattering approach applies not only to conductance 
but also to noise and higher order current correlations, 
the role of inelastic scattering^'* and dephasing*^ in noise 
can also be investigated with this approach. 

A brief account of our main results has been given in 
Ref. 0- We demonstrate that the generating function 
for an arbitrary conductor connected to only one single 
channel voltage probe agrees with the generating func- 
tion of the conductor connected to a dephasing probe. 
A central result of our discussion is that the generat- 
ing function of the one-channel probe models can also 
be obtained by phase averaging the generating function 
of the completely coherent conductor. Phase averaging 
is performed by connecting the conductor with the same 
strength as the voltage or dephasing probe to an exter- 
nal single mode scatterer characterized by a single phase 
with a given distribution'*. Interestingly, the appropri- 
ate phase distribution which leads to a phase averaged 
result identical to that found from the voltage and de- 
phasing probes is simply a uniform distribution. A meso- 
scopic scatterer which provides such a uniform distribu- 




FIG. 1: (color online) A voltage (dephasing) probe described 
by the potential Vp (the occupation function rip), connected 
to a mesoscopic conductor with M terminals. 
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tion is a chaotic cavity. A chaotic cavity, like a voltage 
or dephasing probe, reemits a carrier only with a certain 
time-delay. The importance of the time-delay is mani- 
fest even in the completely incoherent limit: then as we 
demonstratei at least for a particular geometry, the one- 
channel voltage or dephasing probe and phase averaging 
are equivalent to a classical exclusion statistics model^^ 
provided the time-delay is also incorporated into the clas- 
sical model. 

The equivalence of voltage probe and dephasing probe 
models and phase averaging is a special property of one 
channel probes^ Already for probes with more than one 
quantum channels it is known that the voltage probe 
model and the dephasing probe model give different re- 
sults for conductance and noiseJi^^ Furthermore a con- 
ductor connected to two (or more) voltage probes will in 
general exhibit a generating function that differs from the 
one obtained with two dephasing probes.'* We discuss in 
detail why there are differences between phase averaging 
and dephasingiS as soon as there are two or more probes. 

We illustrate our basic results with the help of con- 
ceptually simple conducting structures. A prime exam- 
ple is the Mach-Zehnder interferometer: the simplicity 
here stems from the fact that there are no closed or- 
bits. The electrical Mach-Zehnder interferometer was 
used to discuss de phasing in on e-channel ballistic inter- 
ferometers in Refs. I20l2ll22l23l It has been experimen- 
tally realized in Ref. l24l and has since been the focus of a 
number of theoretical works^SiSSiSLSSiSS followed by ad- 
ditional experiments^Mi. We mention that the recent 
experiment of Litvin et al^ shows a visibility in con- 
ductance as a function of temperature and voltage which 
is in good agreement with theoretical results based on 
voltage probes^. 

Another conceptually simple setup which we will use 
to illustrate our results is a beam splitter structura^SiS^ 
coupled to dephasing and voltage probesi2*iS*^^i2^. This 
example illustrates the dramatic difference between de- 
phasing and voltage probes arising in the case of two 
transport channels: whereas in the presence of a dephas- 
ing probe current-current correlations remain negative, 
a voltage probe can generate positive current-current 
correlationsiLiSi^^, an effect demonstrated in a recent 
experiment by Oberholzer et a\& For a brief discus- 
sion of this and related experiments we refer the reader 
to Ref. m 

In the presence of voltage and dephasing probes deco- 
herence is a consequence of the escape of a carrier out 
of a coherence volumei^. This mechanism differs of 
course from other microscopic dephasing processes, such 
as electron-phonon or electron-electron scattering j22, To 
the extend that voltage (or dephasing) probes are used 
to mimic the effect of these microscopic processes, they 
represent only a phenomenological description. However, 
voltage probes are of course real elements of mesoscopic 
conductors and the escape and reinsertion of carriers are 
true physical processes, if it is taken into account that 
the particles entering the probe are emitted not instan- 



taneously but after a certain delay time. In this respect 
probe models are not phenomenological but provide an in 
principal fully quantum mechanical model which permits 
the investigation of decoherence in a scattering process 
that would in the absence of such probes be completely 
coherent. Indeed several authors have noticed and dis- 
cussed the similarities between the description of inelastic 
processes on the level of quantum kinetic equationa^Siii 
and the physics of voltage and dephasing probes. 

Essentially voltage and dephasing probe models are an 
extension of a coherent scattering problem in which in 
addition to the true current and voltage terminals one 
considers a set of voltage or dephasing probes (see Fig.^). 
Reducing this enlarged scattering probleroi*^ with a set 
of boundary conditions which assure that there is not 
charge accumulation in the probes, one obtains a con- 
duction problem with only the true voltage and current 
contacts. This leads to a description of transport which 
incorporates phase breaking processes to the extend that 
carriers visit the additional probes. On the level of con- 
ductance and noise the reduction from the large coher- 
ent problem to the incoherent fewer contact problem is 
straight forward. The reduction is more challenging if we 
want to capture the entire transport statistics. We for- 
mulate a stochastic path integral approach^i^S*^ for the 
full counting statistics which treats the escape and injec- 
tion process for voltage and dephasing probes. This ap- 
proach uses separation of timescales in an essential man- 
ner: the delay time of the probe has to be much larger 
than the inverse average attempt frequency of scattering 
wave packets. 

The effect of dephasing on the transport statistics for 
entangled states was recently studied, both for a probe 
model with conservation of average current only^ as well 
as for a related model with full current conservation^. 
Attempts at models with dephasing stubs which do not 
generate additional noise has lead to stubs in the form 
of chaotic cavities with long delay timea^. We present 
a detailed discussion of these and other related works in 
section Hxl 



II. VOLTAGE AND DEPHASING PROBES 

Voltage and dephasing probes are used to introduce 
inelastic or elastic incoherent scattering respectively into 
a quantum coherent system. An additional terminal - 
either a voltage or a dephasing probe- is connected to a 
coherent mesoscopic conductor, as shown in Fig. ^ Par- 
ticles entering the probe are later incoherently reemitted 
into the conductor. Scattering via the probe, a parti- 
cle thus looses its phase coherence and in the case of a 
voltage probe, it also changes its energy. 

A voltage probe is a real, physical component used 
in many mesoscopic experimenta^£'^2*^2iSli. It consists 
of a large metallic contact attached to the mesoscopic 
conductor, see Fig. [3 The contact is left floating or is 
connected to a voltmeter, i.e. ideally there is no current 
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FIG. 2: (color online) (a) A voltage probe is a metallic con- 
tact with a floating potential Vp. It can be represented in 
an electrical circuit as a metallic region with capacitance C 
to the ground and charge relaxation resistance R. The low 
frequency current into the probe vanishes, (b) Model of a 
dephasing probe: the probe is described by a non-equilibrium 
occupation function rip. The low frequency current vanishes 
separately in each energy interval. 

drawn at the probe. In response to the injected charge 
on the probe, the potential Vp — Vp{t) of the floating 
probe develops fluctuations on the timescale — RC, 
as sketched in Fig. O Here R is the total charge relax- 
ation resistance from the probe into the M terminals, 
where 1/R — with Gap the conductance from 

the probe to terminal a, and C is the total capacitance 
of the probe. The origin of these fluctuations has a very 
natural explanation: Injected charges raise the potential 
Vp which leads to an increase in the outgoing current. 
This consequently reduces the charge on the probe and 
with this Vp decreases etc. The timescale for the contin- 
uous charging and discharging of the probe is just as in 
classical circuit theory the i?C-time. This picture clari- 
fies, why current and current fluctuations at the probe at 
low frequencies w < I/t^ are completely conserved. Put 
differently, for a measurement during a time much longer 
than Td there is no charge accumulation in the probe. 

We assume that the thermalization in the voltage 
probe is efficient, i.e. charges injected into the probe scat- 
ter inelastically on a timescale much shorter than r^. The 
electron occupation is thus described by an equilibrium 
Fermi-Dirac distribution fp{Vp, T). It is further assumed 
that the temperature T is fixed by the surrounding lat- 
tice, thus the voltage probe is dissipative and does not 
conserve heat current. 

Contrary to the voltage probe, the dephasing probe is a 




FIG. 3: (color online) Time-dependent potential Vp{i) at a 
voltage probe. The current injected from the conductor into 
the probe gives rise to fast fluctuations on the time scale Ts = 
h/eV. The response of the probe gives rise to modulations 
on a much longer timescale Td = RC. 



conceptual tool used to model quasi-elastic dephasing^. 
A particle injected in an energy interval E to E + dE 
is incoherently emitted after spending an average time 
Trf inside the probe, in the same energy interval. This 
allows an energy change of dE much smaller than the 
applied voltage eV or the temperature fc^T. The dis- 
tribution function np{E) in the energy interval is pro- 
portional to the number of carriers in the interval. As 
scattering in each energy interval is independent, np{E) 
is -contrary to the distribution function of the voltage 
probe- a strongly non-equilibrium distribution function. 
Moreover, the fluctuations of the distribution function 
np{E) = np(E,t) in each energy interval are indepen- 
dent. Their origin is qualitatively the same as those of 
the potential fluctuations in the voltage probe: the more 
charges injected into the probe, the larger is np{E). This 
leads to an increase in the outgoing current and thus to a 
reduction of np{E). The fluctuations of np{E) occur on a 
timescale r^, the delay time of the probe. Again there is 
thus no charge accumulation of the probe on timescales 
longer than r^, and the current per energy into the probe 
is conserved up to the frequency l/r^. 



III. FULL COUNTING STATISTICS 

We are interested in the full counting statistics, or the 
distribution of charges transmitted through the meso- 
scopic conductor in Fig. ^ during a measurement time 
T. Let us first consider the conductor without the probe 
when only coherent elastic scattering is present. The 
distribution function is denoted by -P(Q) where the vec- 
tor quantity Q = (Qi, Q2, ■ • ■ , Qm) describes the charge 
transfered into each of the M terminals. P(Q) can be 
expressed in terms of the cumulant generating function 
5(A) by means of a Fourier transformation^ 

F(Q) = j dAe^'^^^-'^-^ (1) 

S{A) = ln^F(Q)e'^-Q. (2) 
Q 

The vector A = (Ai, A2, . . . Am) contains the counting 
variables of the different terminals, the conjugate vari- 
ables to Q. The sum and integrals run over all ele- 
ments of the vector,/ dA = (27r)~*^ / dXi . . . d\M and 
— X]qi Qm- Probability conservation leads to the 
normalization of the generating function S'(O) = 0. 

All irreducible moments are obtained by taking deriva- 
tives of the cumulant generating function with respect to 
the counting variables and evaluated at A = 0. For a long 
measurement time r the transmitted charge into a termi- 
nal a is proportional to r, i.e. one can write Qa = tJq,, 
and the zero frequency cumulants for the current can 
be expressed in terms of the generating function of the 
charge, Eq. Written out explicitely, the first three 
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scattering probabilities. The argument of the logarithm 
in Eq. {Tj) contains the probabilities of all scattering pro- 
cesses, each multiplied by an exponential factor contain- 
ing counting variables that indicate the number and di- 
rection of transfered charges. 

Technically the multi-particle probabilities are ob- 
tained by the expansion of the determinanli^ 



FIG. 4: (color online) A two-particle scattering event from 
channels la and 1/3 into 2/3 and 3/3: the occupation probabil- 
ity is multiplied by the scattering probability and by the expo- 
nential factor containing the counting fields. The scattering 

probability [ det (S^^^'^^^^^ f = |52/3,la'53/3,l/3-'52/3,l/3<S3/5,l<:«P 

represents the two indistinguishable processes shown above in 
the center, with the minus sign originating from the fermionic 
statistics of the particles. 



det[l-^n(A^5^A5-l)] = ^ |det 



cumulants are the average current 

e dS 



{la 



IT dXa ' 

the auto- or cross-correlations 

e2 



and the skewness 



=2 d^S 



d-^S 
i^T dXadXpdX^ 



(3) 



(4) 



(5) 



For a general mesoscopic conductor described by a 
scattering matrix 5, Levitov and LesovikSi derived an ex- 
pression for the generating function 5*0 (where the index 
stands for a coherent system): 



5o 
Ha 



1 

In det 



dt I dE Ha with 



(6) 
(7) 



For a conductor with single mode contacts to the M ter- 
minals, the scattering matrix has the dimensions M x 
M. The matrix h contains the occupation numbers of 
the different terminals with n — diag(ni, n2, . . . , ua/ ), 
and the matrix A introduces the counting fields, A = 
diag(e''^i , e*'^^ , . . . , e*^^^ ) . The generalization to many 
modes is straightforward: the dimension of all matrices 
in Eq. Q grows according to the number of transport 
channels, but importantly all channels in the same ter- 
minal have the same occupation function and counting 
fields. 



Probabilistic interpretation 

In what follows we will frequently use the interpreta- 
tion of the generating function in terms of multi-particle 



{a},{fc} 



fceff)} ie{a} 



Here {a} denotes a set of transport modes -in general 
situated in different terminals-, from which particles are 
injected, and {&} is the set of modes into which parti- 
cles are transmitted. Because of particle conservation, 
the number of elements in set {a} is in each particular 
scattering event equal to the number of elements in set 
{&}. The first sum in Eq. © runs over all possible sets 
{a} and {6}, and represents all possible, distinct ways of 
scattering a number m of particles, with m ranging from 
to the total number of transport modes in the M termi- 
nals. The probability that m particles are scattered from 
{a} to {&} is given by |det |^, where the matrix 

iS|^j^ is formed by taking the intersecting matrix elements 
of the columns corresponding to the elements in {a} and 
the rows corresponding to the elements in {&} from the 
scattering matrix S. The determinant expression origi- 
nates from the indistinguishability of particles and their 
fermionic statistics: all possible transmission processes 
which convert an initial state of m particles in set {a} 
into the final state of m particles in {6} are obtained by 
exchanging successively pairs of the involved fermions. 
An explicit example of a two-particle process is shown in 
Fig.H 

The products over the occupation functions of the dif- 
ferent terminals in Eq. determine the probability that 
exactly the m particles from set {a} are injected. The 
exponent contains the sum of all counting fields of the 
emission channels of set {6}, and the corresponding sum 
of the injection channels of set {a}, multiplied by (—1). 
This exponential factor works as a marker, indicating 
the direction and the number of transfered particles in 
a particular scattering event from {a} to {&} (compare 
also Fig. 

It turns out that in the expression of the expansion 
© independent processes factorize out, while alterna- 
tive processes are added. This explains the integral over 
energy in the generating function Sq — JJ" dt J dE Hq: 
in the long time limit all scattering events at different 
energies are independent, thus the expansion of the de- 
terminant -extended in energy- leads to a product, where 
each factor has the same form but stands for a different 
energy. This can be converted into an integral over the 
function Hq. 
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IV. MACH-ZEHNDER INTERFEROMETER 



The generating function Hq is obtained from Eq. (jS)): 



To illustrate the full counting statistics and the proba- 
bilistic interpretation, we first discuss the Mach-Zehnder 
interferometer (MZI) without a probe. A scheme of the 
electronic MZI is shown in Fig. El It consists of two arms 
connected to four electronic reservoirs 1 to 4 via beam 
splitters A and B. The transmission (reflection) proba- 
bility of the beam splitters is Ta and Tb {Ra and Rb) 
respectively. Transport in the single mode arms is uni- 
directional, corresponding to transport along edge-states 
in the quantum Hall regime. Such a setup was recently 
realized experimentalljiSiiSSiSi. Interference occurs, be- 
cause the electrons have two alternative paths to prop- 
agate through the interferometer between beam splitter 
A and B. An Aharonov-Bohm flux ^ab threads the two 
arms, and the different vector potentials in the two arms 
lead to a phase difference $ ~ ^^^ab- This difference 
creates a characteristic flux-periodicity in the interference 
pattern, the Aharonov-Bohm effect. 



The conceptual simplicity of the interferometer origi- 
nates from the exclusion of backscattering which discards 
closed orbits. As a consequence only the elements in the 
off-diagonal 2x2 blocks of the total 4x4 scattering ma- 
trix S are non-vanishing^*^ . For our purposes we only 
need the processes from terminals 1 and 2 towards ter- 
minals 3 and 4 given in the lower off-diagonal block 



S31 

^32 
^41 



-RaRrc 
TATBe'" 



+ TaTb 
RaRb 



i^RaTb+WTaRbc'^ 
iy/RATBe"^ +i^TaRb. 



(9) 



We consider equal length of the arms, giving energy- 
independent scattering amplitudes. Constant scattering 
phases can be absorbed into The second off-diagonal 
block describes independent processes from terminals 3 
and 4 to 1 and 2 which are not indicated in Fig. [S] and 
which we do not need to specify here. 




FIG. 5: (color online) Mach-Zehnder interferometer: a four- 
terminal conductor with unidirectional transport channels 
(indicated by arrows), threaded by a magnetic flux ^as. 



i/o = ln[((l-ni)(l-n2) 

+ni(l - n2)(T3ie'(^^-^i) + Tiie'^^^-^'^) 
+712(1 - ni)(T32e'(^^-^^) + T42e^(^^-^^)) 
+ nin2e''-^'+^'-^''-^'^)], (10) 



where T^^ — {Sapl'^- Using the probabilistic interpre- 
tation, we can identify the different contributions: The 
term in the first line stands for the case that no particles 
are injected from terminals 1 or 2. The second and third 
line describe one-particle scattering from terminal 1 or 2 
into either 3 or 4. The transmission probability of each 
process is multiplied by the counting factor, for example 
T3ie^^^^~^^^ marks a transfer from 1 to 3. The forth line 
represents two inf ailing particles from 1 and 2, which end 
up in 3 and 4 via an exchange process: the events 1 — > 3, 
2 — > 4 and 1 — > 4, 2 — + 3 are indistinguishable. Its prob- 

abihty is given by | det (s^lll]^ P = \S31S42 ~ 53254ip 
which is here equal to one as a consequence of the Pauli 
principle. 

In the simplest case of zero temperature, and a voltage 
applied at terminal 1, the electron occupation function 
is unity in reservoir 1 and zero in reservoirs 2 — 4 in the 
energy interval < E < eV. Then all charges are inci- 
dent from contact 1, and only one-particle processes take 
place. Each infalling particle has exactly two final states, 
characterized by transmission into either contact 3 or 4. 
Therefore, the interferometer acts like a beam splitter 
and the statistics of transmitted charge is a binomial dis- 
tribution. The cumulant generating function is obtained 
from Eq. ifTUI) as 



^0 = iVln[r3ie 



j(A4-Ai)i 



(11) 



where N = can be understood as the number of 

h 

infalling wave packets during the measurement time r. 
Probability conservation requires T31 -1-741 = 1- Because 
of charge conservation, the generating function depends 
only on the differences A3 — Ai and A4 — Ai. In order to 
obtain all information about the transmitted charge it is 
sufficient to count the number of particles in one of the 
contacts 3 or 4. Here, we set Ai = A4 = and perform 
the Fourier transform with respect to A3. We obtain the 
distribution function PiQ^) = {Qt^x{1 - T^i)'^-'^^ 
for the charge Q3 transmitted into contact 3, a bino- 
mial distribution. Note that the transmission probability 
T31 = RaRb + TaTb - 2^RaRbTaTb cos $ depends on 
the magnetic flux, and the correlations between the cur- 
rent in 3 and 4, C34 = — ^-^T3i(l — T^i), are negative as 
is generally the case for mesoscopic structures in a zero 
impedance external circuit^. 
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Probe coupled to the interferometer 

We want to address the question of how the distribu- 
tion of transmitted charges is affected by couphng the 
MZI to a voltage or dephasing probe? After adding the 
additional terminal as shown in Fig. the setup is de- 
scribed by a 5 X 5 scattering matrix. In the case above of 
zero temperature and a voltage applied only to terminal 
1, the occupation functions of terminals 1 — 4 are known 
but the occupation of the probe Up is not determined. 
While in the coherent interferometer only one-particle 
processes with two possible final states were allowed, the 
additional probe opens up additional paths and leads to 
two-particle processes. It is instructive to group the dif- 
ferent contributions to Hq into scattering processes ac- 
cording to the number of particles entering or leaving 
the probe. This gives 



Hn 



In [(1 - Up) {{RaRb{1 ~e) + TATB)e'^' 
+ {RATB{l-e)+TARB)e'^' 
+2^RaTaRbTb{1 - £)cos$(e^^^ - e'^^) 

+np {{RaRb + TaTb{1 - £))e*^^' 
+ {RATB+TARB{l~e))e 



+2y/ RaTaRbTbII-e) cos 'l>{e 



i\4 



+ npeTAe''-^-'+^'-^''^ -t- (1 - np)sRAe' 



,(12) 



The different processes are illustrated in Fig. |(7|) and con- 
sist of i) no charge moving into or out of the probe (the 
first three lines in Eq. p2fl . multiplied by (1 — rip)), ii) 
processes where one particle is entering and one is leaving 
the probe (line 4 to 6, multiplied by rip), iii) processes 
where one particle is leaving the probe (the first term in 
the last line) and iv) one particle is entering it (the very 
last term in Eq. ^^). 

As discussed in section ^1 a defining property of the 
probe is charge conservation on timescales longer than 
the delay time Td of the probe. This means that all cu- 
mulants of the current flowing into the probe are zero 
for low frequencies uj ^ l/rd- Clearly, the expression for 




FIG. 6: (color online) The upper arm of the Mach-Zehnder 
interferometer is coupled to a dephasing or voltage probe with 
coupling strength e. 



the generating function, Eq. H12() is then not the solution, 
since current and current fluctuations at the dephasing 
probe obtained from derivatives with respect to Xp are 
non-zero. Also the occupation Up (and Xp) is fluctuating, 
which is not taken into account in Eq. (|12() . Therefore, 
this equation is defining an unconstraint generating func- 
tion: it treats the dephasing probe like any other charge 
absorbing terminal. The question is now, what is the cor- 
rect expression for the FCS for a dephasing or a voltage 
probe attached to the conductor? 



V. STOCHASTIC PATH INTEGRAL 

To arrive at the correct answer we make use of the 
stochastic path integral (SPI) technique developed in 
Refs.|42,J3f The stochastic path integral approach is one 
of a number of technics which can be used to find the 
generating function&S. The cornerstone of the stochas- 
tic path integral technique is the separation of timescales. 
In mesoscopic systems the current typically exhibits fluc- 
tuations with two qualitatively different physical origins: 
i) the intrinsic current fluctuations, resulting from prob- 
abilistic scattering of individual charges, which occur on 
a timescale of the order Tg, and ii) the fluctuations of 
potentials or distribution functions arising as a current 
conserving response to the intrinsic fluctuations, which 
occur on a timescale of the order r^, the dwell or delay 
time (compare also Fig.Ol. 

For our system, a general mesoscopic scatterer con- 
nected to a voltage or a dephasing probe (see Fig. ^1, 
the short timescale is given by the average distance in 
time between two successive wave packets, typically of 
the order of — h/eV. The long timescale is given by 
the delay time of the probe r^. It is thus essential for the 
application of the SPI that 



Ts < Td 



(13) 



i) U 





FIG. 7: (color online) i) No particles are moving in or out 
of the probe, ii) One particle enters, one leaves the probe, 
thus there is no net charge transfer into the probe, iii) One 
particle is emitted from the probe, iv) One particle enters the 
probe, i) and iv) are one-particle processes, and ii) and iii) 
two-particle processes. Note full red dots represent injected 
particles and filled (empty) green circles represent the one 
(two) final states. 
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or in other words that the apphed voltage is large enough 
eV > h/Td. 

Based on the separation of timescales, we introduce 
for the derivation of the SPI an intermediate time step 
Ai with Ts ^ At ^ Td, during which the charge in the 
probe (and its counting variable Ap) changes marginally. 
The FCS for this time interval is given by the Levitov- 
Lesovik formula Eq. © for a scatterer connected to the 
probe, i.e. the unconstraint generating function. The 
charge fluctuations in two successive time intervals Ai 
can be expressed in terms of the distribution functions 
of charge transfered during each of the two time inter- 
vals. By stepwise extending this result to longer times 
and taking the continuum limit, we finally obtain the 
distribution function -P(Q) in terms of a stochastic path 
integral. A formal derivation of the SPI is given in the 
appendix 1X1 here we simply state the results, first for the 
dephasing probe and thereafter for the voltage probe. 



A. Dephasing probe 

In the dephasing probe, the low frequency currents 
per energy interval into the probe are conserved due to 
fluctuations of the occupation number np{E,t). For a 
coherent conductor to which a dephasing probe is con- 
nected, the generating function S in Eq. ||2Jl is obtained 
by the stochastic path integral over all possible distribu- 
tion functions Up and values of the counting variable Ap 
at the dephasing probe, 



exp(S') = / 'DnpVXpexp{S), 



where 



^4 



dt / dE [—iTdXpTip + Hq] 



(14) 



(15) 



Thus the additional term —iTdXpiip in Eq. H15|) takes into 
account charge conservation on the probe. 

In general the path integral in Up and Ap can not be 
evaluated exactly. It is however possible to evaluate the 
integral in saddle point approximation. The underlying 
idea of the saddle point approximation is that the main 
contribution to the integral comes from the region around 
the stationary, or saddle point, where 



6S 
5n„ 



SS 

SXr, 



(16) 



Note that here both rip and Ap depend on both energy 
and time. Away from the saddle point the integrand 
fluctuates rapidly as a function of Up and Ap, giving a 
small contribution to the integral. We thus expand S to 
second order around the saddle point as 



S 



1 A .2^^^ 



1a 2-^'^ 



AX A 

'AXpAup 

oAoOno 



(17) 



Here AAp = Ap — A^ and Aup 



Up are solutions to the saddle point equations H16(l . The 
zeroth order solution S° is just S with the saddle point 
solutions Ap and rip inserted. As is shown in the appendix 
El the Gaussian fluctuations around the saddle point only 
give rise to corrections of the order Tg/rd <C 1. This 
means that as long as the underlying condition for the 
SPI, Eq. H13|l is fulfilled, the saddle point approximation 
S" = S"" is a good one. 

Carrying out the functional derivatives 6S/Snp and 
6S/5Xp = we arrive at the saddle point equations for 
"p 



Ap and 7i. 



1 dHo 
Td dXp 



iXr, 



1 dHo 

Td dup 



(18) 



In this paper we will only consider the stationary limit, 
T ^ Td- In this limit the time derivatives Ap and fip 
can be neglected, the functional integrals are reduced to 
normal integrals and the function S = S becomes pro- 
portional to the measurement time 



= - / dEHo 



(19) 



Consequently, the saddle point equations are reduced to 



dX„ 



= 0, 



dHo 
drir, 



= 0. 



(20) 



In the general case, with many transport modes in the 
contacts between the probe and the conductor, the saddle 
point equations are highly nonlinear in both Ap and rip 
and it is only possible to solve them numerically. In the 
case with a single mode the solution can be found analyt- 
ically. Just as for the special case of the MZI discussed 
above (see Fig.|7|l, all possible scattering processes belong 
to one of the four different groups: i) no particle injected 
into and no particle emitted from the probe, ii) one par- 
ticle injected into and one emitted from the probe, iii) 
no particle injected but one emitted and iv) one particle 
injected into but none emitted. Due to the Pauli prin- 
ciple no other processes exist, i.e. it is not possible for 
more than one particle to leave or enter the single mode 
probe in the same scattering process. Making use of the 
expansion of the generating function in terms of multi- 
particle scattering probabilities, Eq. we can write the 
generating function as 



Ho{np, Ap) = In [(1 - Up) (goo + goie'^") 



gii + gioe 



~iXr 



)] 



(21) 



Here we introduced the notation qa for the total multiple 
scattering probabilities (multiplied with the appropriate 
counting field expressions). The index Z = 0, 1 denotes 
the number of particles injected into the probe, and k = 
0, 1 the number of particles emitted from the probe. It 
is possible to express q^i in a compact form in terms of a 
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number of determinants F^i, defined as 

Fki = det[l + nfc(Al5U/5-l)] 
Uk = diag(ni,?i2, ...nM,/c) 



Afc = diag(e^^\e^ 



,k). 



(22) 



Here, and Xk are diagonal matrices of the distribution 
functions and counting fields respectively, with the last 
entry, corresponding to the probe terminal, just given by 
the number fc = 0, 1. The relation between qki and Fki is 
given by 



900 — ^00 7 901 — -foi — Fqo, 
qio = Fio, qii = Fii — Fio- 



For Hq in Eq. (|21|l . the second saddle point equation in 
(Pn|l gives 



dHo 



= - (goo + 9016'^") + qn + qwe = (24) 



from which one obtains the expression for A!J 



qoo 



^(gii - 9oo)^ + 4goi'7io 



2gi 



01 



(25) 



The correct sign in this solution is defined by the nor- 
malization Ho{A = 0) = 0. Inserting this back into the 
generating function, noting the terms proportional to Up 
drop out, we get the generating function 



S ^ Nln 



46 



(26) 



with A = qoo + qii = Fqo + Fu — Fm and B = 
qmqii - goigio = Fo^Fh - FoiFio and N = eVr/h. 
This is the general expression of the constraint gener- 
ating function for a conductor coupled to a single mode 
dephasing probe. The interpretation of this function is 
not as straightforward as the probabilistic expansion of 
the unconstraint function Hq, Eqs. {Tj) and ||SJ). Both 
the parameters A and B contain multi-particle scatter- 
ing processes, but A represents only processes without 
any net charge transfer into the probe. The square root 
term technically originates from the fact that the sad- 
dle point equation Eq. is quadratic in e*"^", which 
reflects that in one scattering event, the charge in the 
probe can be increased or diminished by one. 



B. Voltage probe 

In the voltage probe, the potential Vp fluctuates in 
response to the injected current fluctuations. Just as 
for the dephasing probe the potential fluctuations lead 
to current conservation on the timescale of t^. For the 
voltage probe it is however the total energy integrated 
current that is conserved. The effect of the fluctuating 
voltage on the FCS can again be incorporated via the 



stochastic path integral over the potential Vp and the 
energy independent counting variable A„ as 



exp(S'y) = / VVpVXpexj>{Sv), 



where 



dt 



ir^XpcVp 4- / dEHi 



(27) 



(28) 



In the long measurement time limit t ^ Td the time 
derivatives can be neglected and we have, analogous to 



(23) Eq. CH) 



dEHn, 



(29) 



where we note that Vp enters in the distribution function 
Up of the probe. The saddle point equations are 



dXn 



0, 



dSy 

dVr, 



(30) 



Focusing on the zero temperature limit and on energy 
independent scattering only, we can write the function 



er 



eVp i-eV 

Tip — 

JeV„ 



{VpHo\n, = l + iV- Vp)Ho\n,=o) ■ (31) 



This function is linear in the applied voltage V and in 
the potential of the probe Vp even for many transport 
channels. The second saddle point equation in Eq. 
makes the part proportional to Vp drop out and at the 
same time determines the saddle point solution for the 
counting field A^. For the constraint generating function 
we obtain* 



S\, 



cVt 



= 1,A0- 



(32) 



For a single transport mode the solution to the second 
saddle point equation, Eq. (|30|) , is identical to the one for 
the dephasing probe and we arrive again at the expression 
in Eq. The generating function Sy is consequently 

identical to Eq. H26|l , Sy = S. There is thus no difference 
between a dephasing and a voltage probe for the case of 
a single mode probe and energy independent scattering. 
However, if there are more than one single mode probe, 
multi-mode probes or energy dependent scattering, the 
different electron occupation in the voltage and dephas- 
ing probes lead to different transport statistics. This is 
explained in section IVIIII for an example of a conductor 
with two transport channels coupled to either a voltage 
or a dephasing probe. 
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VI. THE MACH-ZEHNDER 
INTERFEROMETER WITH THE PROBE 



The coherent Mach-Zehnder interferometer at zero 
temperature and with a voltage V apphed at terminal 1 
is characterized by a binomial distribution, as explained 
in detail in section HVI We first note that the cumulant 
generating function 5*0, Eq. Hll|l . can be rewritten as 



Sq = iVln[6 + ccos<I>] 



(33) 



Here, the different terms in the transmission probabil- 
ities T31 and T41 are rearranged into the parameter b 
and c. The parameter b = {T/Tb + RaRb)^'^'' + 
{TaRb + RaTb)^^^^ represents the classical contribu- 
tion due to particles which go either along the up- 
per or the lower arm, and ccos$ stands for the co- 
herent quantum interference contribution with c = 



2^/RaTaRbTb (e^ 
is set to zero. 



The counting variable Ai 



Fig. shows the MZI with a dephasing probe con- 
nected to the upper arm. The full counting statistics for 
the interferometer with the probe attached is obtained by 
applying the formalism of section FV AI Solving the sad- 
dle point equations in Eq. (^01) with Hq given by Eq. ifT^ 
we find for the constraint cumulant generating function 



iVln 



Ml 



e cos $ 



(34) 

In comparison to the result for the coherent interferom- 
eter, Eq. (|33(l . the interfering contribution is attenuated 
by the factor \/l — Thus, as expected the generating 
function becomes flux independent in the limit of strong 
dephasing 



case of equal beam splitters, Ra = Rb = R = I — T 



C33 — 
C333 = 



-{R^ - 2i?TVT^cos$), 
— '2RT [{R^ - RT + T^) 



+ {R - r)2vT^cos$ - RT{1 - e) cos 2$) 
— '^RT [-{R - Tf{R^ - 3i?T + T^) 



(36) 

(37) 

)]- 
(38) 



+{2R^T'^-{R-Tf{R^-mT+T'^))Vl^cos^ 
+3RT{R - Tf{l - e) cos 2$ 

-2i?2T2(l--e)icos3$ . 

Each cumulant contains oscillating terms proportional to 
cos($), cos(2$), . . . ,cos(fc$), where k is the order of the 
cumulant. As is clear from Eqs. H36|) - (|38|l . the effect of 
coupling the MZI to the probe is to multiply the oscillat- 
ing terms proportional to cos(fc$) with damping factors^^ 
(1 — t)^/'^. Importantly, the ^-independent terms in the 
cumulants are however not affected by the coupling to 
the probe. Despite the complicated multi-particle origin 
of the contributions to the generating function, Eq. JSiJ, 
the only effect of the probe is thus to damp the AB- 
oscillations in the cumulants. 

Similarly, a second single mode probe attached to the 
MZI with coupling strength e' results only in modifying 
the damping factor of the oscillations in the cumulants 
in the following way: (1 - e)''/^ ^ (1 - £)''^^(1 - e')''^^- 
Again it does not matter if the additional probe is cou- 
pled to the same or the other interferometer arm. 

For energy-independent scattering considered here, the 
discussion at the end of section IV Bl applies and Eq. H34() 
does not change if the dephasing probe is exchanged with 
a single mode voltage probe. This is however not the case 
for more than one probed. 



VII. PHASE AVERAGING 



Se=i = N\n 



(35) 



In addition, the parameter b is in Eq. H34(l multiplied 
by a reduction factor (1 — £/2), and a third contribu- 
tion ev6237i2/2 appears. As pointed out in section Hvl 
the processes into and out of the probe give rise to ex- 
change, or two-particle interference processes, for exam- 
ple 1 — !■ 3,p — > p and 1 — > p,p — > 3, which are indistin- 
guishable. Note that the processes with particles moving 
along either of the arms (6) are reduced, regardless of the 
fact that the probe is coupled to the upper arm only. 

The first two cumulants obtained from Eq. H34|) re- 
produce known results for current and nois822*2^iS&. For 
completeness we present the first three cumulants for the 



The reduction of the oscillating terms in the cumulants 
due to the coupling of the MZI to a probe, indicates that 
the effect of a dephasing probe can be explained in terms 
of phase averagingSSiS^. Consider the interferometer cou- 
pled with strength e not to a dephasing probe, but to 
an additional elastic coherent scatterer (p as depicted in 
Fig. |S1 Every particle incident on contact (p is returned 
coherently with a phase factor e"^. An electron in the 
upper arm of the interferometer can enter the contact 
<j) and after multiple internal reflections continue on its 
path with an additional phase factor e^"^ . This additional 
phase factor enters into the coherent generating function 
as ifo(*) ->i/o(* + '^)- 

We make the simplest possible assumption that the 
phase ^p is uniformly distributed between and 27r. Then 
the electron phase (p is related to the scattering phase by 
(p = (p + -K + 2 arctan[Vl — esm (p/{l — Vl — £ cos (p)] and 
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obeys a periodic distribution function^ 



1 dip 
27r dip 



1 



2tt2- 



■ + 2\/l — ecos{(p) 



(39) 



The average of the coherent generating function for the 
MZI is with Eq. ^ given by 



2tt 



(40) 



and agrees exactly with the result (|34|l obtained from the 
dephasing probe model. 

Interestingly, it is possible to obtain such a phase av- 
erage of the cumulant generating function by taking the 
additional elastic, coherent scatterer to be a single mode 
chaotic cavity with a long dwell time, equal to the delay 
time of the probe ^ h/eV. The scattering phase ip 
picked up by electrons scattering at the cavity depends 
on energy. It is known from random matrix theorjsi^ 
that scattering amplitudes at two different energies E 
and E+dE are uncorrelated if the energy difference dE is 
much larger than the inverse dwell time h/rd- Put differ- 
ently, en electron at energy E sees a completely different 
cavity scattering potential than an electron at i? -I- dE. 
Since eV 3> h/rd, the integral over energy of the cumu- 
lant generating function, Eq. (|19|l . for the MZI connected 
to the cavity effectively means a sum over a large num- 
ber of realizations of cavity scattering potentials, i.e. an 
ensemble average. It is also well known from random ma- 
trix theorjii^ that the phase ip of an ensemble of single 
mode cavities is uniformly distributed between and 27r. 
Formally, since the scattering in the MZI itself is inde- 
pendent on energy in the interval to eV, we can thus 
write Jq dE 



lo^ '^'^^ ^'^"^ arrive at the phase 



average of Eq. 

We note that from Eq. H40|l it is clear why the de- 
phasing probe model exclusively affects the oscillating 
contributions to the cumulants: the function _ffo can be 




FIG. 8: (color online) Phase averaging: a particle enters the 
elastic scatterer with probability e and obtains a uniformly 
distributed phase factor e"^. After multiple internal reflec- 
tions it leaves the scatterer with an overall phase e"^. The 
inset shows a general conductor to which an elastic scatterer 
A is connected. 



expanded in terms of the cumulants, and the phase av- 
eraging leaves the contributions that are independent of 
<f> unchanged while it damps the oscillations. 



General conductor 

The correspondence between phase averaging and the 
dephasing/voltage probe model also holds for a general 
mesoscopic conductor with M terminals. Instead of a 
single mode (dephasing/voltage) probe, an elastic scat- 
terer (p is via a single mode connection coupled to the 
conductor as shown in the inset of Fig. |S1 

We define a (M -I- 1) x (M -I- 1) scattering matrix U of 
the conductor with the additional contact (j) that can be 
written in a block form 



U 



(41) 



Here the AI x M block Sqq describes scattering between 
the M terminals of the conductor, the M x 1 (1 x M) 
block 5^0 (Sq^,) the scattering from (to) the conductor 
to (from) the phase contact (j) and complex number 
the reflection from the contact (p back to (p. The total 
M X M scattering matrix S for the conductor connected 
to the elastic scatterer 6 can then be written as 



S = S, 



00 



00 



1 - S^e'^ 



■00- 



(42) 



Inserting this scattering matrix into Eq. the generat- 
ing function Ho{(p) depends on ip. We then average the 
cumulant generating function Ho{ip) over the uniformly 
distributed phase (p as 



(43) 



We show in appendix ^ that the averaged generating 
fmiction (Hq)^ is identical to the result for the dephas- 
ing (or voltage) probe, Eq. (|26(l . valid for an arbitrary 
mesoscopic conductor connected to a single mode probe. 

It is important to point out that for any geometry con- 
taining two or more single mode probes, the phase aver- 
aging approach does in general not give the same result 
as the dephasing probe model. The reason for this is 
that interference between scattering paths that pass the 
different probes/contacts Pj/4>j in different order (but 
the same number of times) survives phase averaging but 
is suppressed by dephasing probes. A more detailed dis- 
cussion of this together with concrete examples where the 
dephasing probe and phase averaging approaches give dif- 
ferent results, are presented in appendix^l We also point 
out that for probes with more than one mode there is no 
obvious generalization of the phase averaging approach. 
An investigation in various models for phase averaging 
in multi-mode probes is beyond the scope of the present 
paper and thus deferred to a future publication. 
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VIII. TWO-CHANNEL STRUCTURE 

An interesting system where the couphng to voltage 
and dephasing probes strongly affect the transport prop- 
erties is shown in Fig. |3 This system was investigated 
theoretically by Texier and BiittikeriiLi^ and recently re- 
alized experimentally by Oberholzer et al& It consists 
of a three terminal quantum hall bar to which an addi- 
tional probe is connected. The coherent setup without 
the probe is displayed in the inset of Fig. |5| The filling 
factor is tuned to two, such that the current is carried by 
two parallel edge states propagating along the boundary 
of the sample. 

The two quantum point contacts QPC 1 and QPC 3 
work as beam splitters with Ti (Ri) and T3 {R3) the 
transmission (reflection) probabilities respectively. Both 
quantum point contacts are completely open for the outer 
channel. We consider the case of zero temperature and 
a voltage is applied at terminal 1, while the terminals 2 
and 3 are kept at ground. 

In the original investigation the focus was on the cross- 
correlations between the currents in terminal 2 and 3. 
It was shown by Texier and Biittiker that the sign of 
the cross-correlations depends on whether a voltage or a 
dephasing probe is coupled to the system. In the case 
= 0, the inner channel leads to contact 2 and the 
outer to contact 3. Then for the coherent system the 
cross-correlations C23 vanish, a voltage probe introduces 
positive cross-correlations, and for a dephasing probe the 
correlations turn out to be negative. The positive cross- 
correlations have been measured recently by Oberholzer 
et al.^. 

It is interesting to understand these statistical effects 
on the level of the generating function. The generating 
function especially enlightens the difference between volt- 
age and dephasing probes, which is present only in the 
case of more than one transport channel if scattering is 
energy independent. 

We start again with the coherent system without the 
probe. From Eq. © the generating function Sq reads: 

50 = A^ln[e*(^^"^i) (44) 
x(l -Ti +rir3e'(^-^-^i) +Tii?3e*(^^-^i))] 

The factor e*''''-'^'*'!^ stands for the outer channel, where 
particles are transmitted from 1 to 3 with probability 
one, it represents a noiseless contribution to the current 
into terminal 3. The factor in the parentheses represents 
the inner channel, where an infalling particle is either 
reflected at beam splitter QPC 1 or transmitted into ter- 
minal 2 with probability T1R3 or into 3 with probability 

Tir3. 

For T3 — 0, the inner and outer channel lead into 
different contacts 2 and 3. Since in the absence of in- 
terchannel scattering the two modes are uncorrelated, 
the cross-correlations vanish, C23 = 0. When T3 is 
greater than zero, the cross-correlations are negative, 
C23 = — ^-^T'^T^R-i, because QPC 3 acts as a beam split- 



ter for the inner mode (but the inner and outer channel 
are still uncorrelated). 

We then consider an additional probe perfectly coupled 
to the system, see Fig. |^ To keep the discussion simple, 
we concentrate on the case T3 = 0. The unconstraint 
function Hq can be written in terms of independent one- 
particle scattering processes 



He, = In 



gi(Ap-Ai)^2 — Ti + Tj^e^^-^p^-^i)) 



(45) 



X (1 - np + npe^(^^"^''))(l - Up + npe'^^^"^"'^) 



The first two factors correspond to the two transport 
modes emerging from terminal 1, and the rest describes 
the inner and outer channel from p to 2 or 3 respectively. 



A. Voltage probe 

In the presence of a voltage probe, rip is an equilib- 
rium Fermi function, and the potential Vp fluctuates in 
order to maintain zero net current flow into the probe. 
The formalism to obtain the FCS is presented in section 
IV Bl Because the current is transported by two modes, 
the argument of the logarithm in Eq. H45|l is quadratic in 
Up, but the energy integrated function Sy at zero tem- 
perature, Eq. H31|l is still linear in Vp. The saddle point 
equations, Eq. H3UII are easily solved, and the cumulant 
generating function is 



Nhi[e- 



A2+A3 



(l + Ti(e' 



A2+A3 



(46) 




FIG. 9: (color online) Beam splitter with a probe: a three 
terminal quantum hall bar at filling factor 2 and one volt- 
age/dephasing probe. The two quantum point contacts QPC 
1 and QPC 3 are open for the outer transport channel. A 
voltage is applied on terminal 1, thus the channels shaded in 
(light) grey are effectively empty and do not matter. The in- 
set shows the setup for the case when the connection to the 
probe is closed, i.e. when there is no probe connected to to 
the beam splitter. 
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The first factor represents a noiseless stream from 1 to 
p which is divided into two noiseless streams in the two 
outgoing channels from p leading into contacts 2 and 3. 
The part inside the parentheses describes the inner chan- 
nel. Passing through QPC 1 the incident stream at the 
probe is noisy and it is again equally divided into the 
outgoing channels. 

Both outgoing channels are fully filled up to the energy 
eVp = eVp{t). This creates a perfect correlation between 
the currents in 2 and 3, hit) = hit), and is the reason for 
the positive cross-correlations given by C23 = ^7r""^^T^- 
Technically the positive sign can be seen in the generat- 
ing function Eq. (|46|l . where the counting fields A2 and 
A3 appear only as a sum in the exponent. For a finite 
transmission T3 the inner edge is split at the QPC 3, and 
the negative cross-correlations of this splitting competeii 
with the positive correlations induced by the potential 
fluctuations at contact p. 

The voltage probe thus acts as a noise divider^, be- 
cause both currents in 2 and 3 show fluctuations, in con- 
trast to the coherent system where only the inner channel 
leading to contact 2 was noisy. As a consequence of cur- 
rent conservation, the noise of the current incident at 
the probe -created by the beam splitter QPC 1- equals 
the sum of current correlations emerging from the probe: 
^RiTi = C22 + C33 + 2C23. Note that both auto- and 
cross-correlations vanish for Ti = separately. 

The counting fields A2 and A3 appear with a factor 1/2 
in Eq. We emphasize that the factor 1/2 does not 

mean that the voltage probe emits half particles with 
charge e/2. The generating function only makes sense 
in the long time limit when many charges are emitted. 
Therefore, the factor 1/2 simply describes the equal par- 
tition of the incident charges into two outgoing channels. 



B. Dephasing probe 

If instead of a voltage probe, a dephasing probe is cou- 
pled to the system, the cross-correlations are manifestly 
negative. The argument of the function Hq, Eq. is 
quadratic in the occupation number Up, therefore, the 
saddle point equations, Eq. (|20|l are nonlinear. For the 
case Ti = 0, when only the outer fully filled channel is 
transported to the probe, the saddle point equations are 
solved by = 1/2 and A° (A2 + A3)/2. The resulting 
generating function is 

5T,=o = iVln[i(e^^+e'^)2] (47) 

Again the counting fields A2 and A3 appear with a fac- 
tor 1/2 displaying equal partition of the incident charges 
into the outgoing channels. In contrast to the voltage 
probe, the occupation of the channels going out from the 
dephasing probe is non-unity. This leads to noise in each 
of the two outgoing channels, despite the fact that the in- 
coming current for Ti = is noiseless. As a consequence 



of current conservation = C22 + C'33 + 2C23 the cross- 
correlations are negative: with C23 = = — C22 — 

~C'33- 

We have not been able to find analytically the exact 
constraint cumulant generating function for arbitrary Ti . 
By expansion in the counting fields A2 and A3 the saddle 
point equations can be solved separately for each order in 
the counting fields and the solutions are found cumulant 
by cumulant. Another method is the expansion in small 
Ti which leads to the constraint generating function to 
second order in Ti 

+ 2(2 + ri)-4e^^^(l + ri)) +0(Tf) (48) 

The cumulants up to second order obtained by this gen- 
erating functions are exact, because they are of second 
order in Ti — 1 — The cross-correlations are al- 

ways negative and given by C23 = and the 

auto-correlations are C22 = C33 = ^T~^('^ ~ 3-Ri). 
The noise incident on the probe is ^-^Ri{l — i?i) — 
C22 + C33 + 2C23. Note that the incident noise is max- 
imal for Ri — 1/2, while the noise in terminal 2 and 3 
has its maximum at i?i = 3/4, and the cross-correlations 
are monotonously decreasing with Rf. 

IX. COMPARISON TO OTHER MODELS 

It is interesting to compare our approach to closely re- 
lated work. An important subject is phenomenological 
models of dephasing in chaotic cavities. Polianski and 
Brouwer— investigated general time-dependent scatter- 
ing in a chaotic dot, using a model where the dot is 
connected to a short stub with a fluctuating potential. 
Although the fluctuating potential leads to dephasing, it 
also creates noise which is an undesired feature for model- 
ing pure dephasing. To overcome this problem Beenakker 
and Michaelis recently developed a stub model''^ with a 
long stub in which the dynamics was chaotic, i.e. in fact 
representing the stub as a second chaotic cavity. In the 
limit of a long dwell time of the stub they found that 
the additional noise created by the fluctuating potential 
in the stub was negligible compared to the regular shot 
noise. 

Beenakker and Michaelis compared their model to a 
dephasing probe model with a static distribution hmc- 
tion of the probe, chosen to yield zero average current 
into the probe. In such a dephasing probe model, the 
low frequency probe current in each measurement is not 
zero. Importantly, in our work we instead consider a 
probe model where the distribution function fluctuates 
in time as a response to the injected charge. The low 
frequency current is thus conserved in each measurement 
and intrinsically no additional noise is created. Conse- 
quently, our probe model has both the desired features of 
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a pure dephasing model, it does not introduce any addi- 
tional noise and it conserves charge, or equivalently the 
low frequency current. It is at the present not clear to us 
whether an extension of the model of Ref. 46 to the FCS 
would yield the same result as our dephasing probe model 
in the corresponding parameter (many modes, weak cou- 
pling) limit. 

In another recent work^ San- Jose and Prada devel- 
oped a method for investigating the effects of a charge 
conserving probe on the full counting statistics of co- 
herent mesoscopic conductors with energy independent 
scattering. They consider a set of localized wave packets 
-"scattering events"- incident on the conductor during 
a time step At = h/eV. In the absence of probes, scat- 
tering events at different time steps are independent, and 
the Levitov-Lesovik result, our Eq. ||HJ), is recovered. In- 
troducing probes connected to the conductor, scattering 
at different time steps becomes correlated: the number 
of charges injected into the probe at one time step de- 
termines the number of charges emitted from the probe 
in the next step (or next few steps). The full count- 
ing statistics is obtained in a way qualitatively similar 
to exclusion statistics models^. In contrast to the stan- 
dard exclusion models the probabilities in the scheme of 
Ref. are multi-particle scattering probabilities. 

In the model of Ref. all processes happen within the 
timescale At. The excess (deficit) charge in the probe 
is at most a few electrons. There is no separation of 
time-scales. This is clearly different in construction from 
our model where the probes contain a large number of 
electronic states and consequently respond to injected 
charges only on the timescale Td ^ At = Ts- Never- 
theless, in the long time limit, we find that the FCS for 
an arbitrary mesoscopic conductor connected to a single 
mode probe is given by our Eq. H26|l . which holds both 
for dephasing and voltage probes. 

In the case with multi-mode probes we have not found 
any general relation between the method of Prada and 
San Jose and our method. It should however be noted 
that for a number of simple examples investigated, we 
find that the method of Ref. and our voltage probe 
(but not the dephasing probe) approach give the same 
long time transport statistics. 

A third related work is by Jakobs, Meden, Schoeller 
and Enss about phase averaging versus dephasing in a 
ballistic, one-dimensional, three barrier setup-- . They 
calculate the conductance as a function of temperature. 
Finite temperature leads to phase averaging due to the 
energy dependence of the phase picked up in scatter- 
ing processes with multiple loops between the barri- 
ers. They find that the phase averaged conductance is 
G ~ VG1G2G3 in the limit of small conductances Gi of 
the individual barriers. 

This is different from the result obtained by summing 
up the classical probabilities for the scattering paths 
through the system. Clearly, this latter result is also 
obtained by coupling each of the inter-barrier regions to 
a dephasing probe. ^ The different results can, along the 



same line as at the end of section IVIII be explained in 
terms of the different effect of phase averaging and de- 
phasing probes on interfering scattering paths that pass 
the two probes in different order. A detailed discussion 
of this is presented in appendix IeI 



X. CONCLUSION 

In this work we have developed a theory of electron sta- 
tistical transport for conductors with voltage and dephas- 
ing probes. Using a stochastic path integral approach we 
have formulated generating functions which determine 
conductance, noise and all higher order current cumu- 
lants. The extension of the discussion of probe models 
for current and noise to the level of full counting statistics 
makes it possible to compare the different models with 
each other on the level of the generating functions. The 
validity of this approach rests on the fact that probes 
respond to a change in their charge state with a certain 
delay time similar to a real voltmeter. Therefore, cur- 
rent conservation at these probes is not instantaneous: 
only the low frequency current cumulants into the probe 
vanish. 

We find that for one single mode probe the generating 
functions for a conductor attached to either a dephasing 
or a voltage probe are identical. They are also equiv- 
alent to the phase average of the coherent generating 
function over an appropriate phase distribution, given 
by the phase distribution of a chaotic cavity. We empha- 
size that these probes are connected to the conductor in 
the same way: it is only the "boundary" condition at the 
probe which differs. These boundary conditions lead to 
an equilibrium Fermi distribution for the voltage probe, 
a non-equilibrium distribution for the dephasing probe 
and leads to coherent reflection in the case of a phase 
probe. In the phase averaging probe used here we av- 
erage over a very specific phase which is external to the 
original coherent conductor. It is important that these 
basic facts are kept in mind when assessing the finding re- 
ported here that all three probes are in the single channel 
limit equivalent. 

For more than one probe, (or a probe with two quan- 
tum channels) however, the equivalence does not hold. 
For the phase averaging procedure used here, certain in- 
terference terms can survive. Thus phase averaging and 
dephasing probes give different results. This is reminis- 
cent to the well known fact that in mesoscopic conductors 
weak localization survives phase averaging but dephasing 
suppresses weak localization. Similarly, phase averaging 
probes leave interferences between certain paths unaf- 
fected but dephasing probes destroy the interference be- 
tween these paths. 

In the case of multiple or multichannel probes it is not 
only the equivalence between phase averaging probes and 
dephasing probes which is lost but also the equivalence 
between dephasing and voltage probes. The different 
electron distributions in dephasing and voltage probes 



14 



become important and lead to different transport statis- 
tics. A spectacular emergence of this is the opposite sign 
of the current cross-correlations caused by either a volt- 
age or a dephasing probe coupled to a conductor with 
two channels. 

The full counting statistics of conductors in the pres- 
ence of dephasing or voltage probes provide a deeper un- 
derstanding of the physics of these probe models and will 
hopefully be a useful point of reference for future inves- 
tigations. Many additional questions are of interest, for 
example how the probes affect the frequency dependent 
cumulants, or the connection between the semiclassical 
stochastic path integral and a full quantum approach or 
simply an investigation of the limit of large numbers of 
probes connected to a mesoscopic conductor. 
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APPENDIX A: DERIVATION OF THE 
STOCHASTIC PATH INTEGRAL 



The FCS for a conductor with M terminals connected 
to a dephasing probe can be derived in terms of a 
stochastic path integral. Since the charge on the de- 
phasing probe is conserved for each energy interval dE 
and the occupation function np{E,t) fluctuates indepen- 
dently in each interval, the distribution of charge trans- 
fered through the whole conductor is independent in each 
energy interval. Therefore, in this section we consider the 
charge transported in a small energy interval dE. 

First, an intermediate time interval At is considered 
with Ts <C At ^ r^. We recall that — eV/h is the 
rate at which carriers injected into the conductor and t^^ 
is the much slower relaxation rate of distribution func- 
tions. The inequality At <C Td implies that the charge 
transfer during time At only results in a marginal change 
of the charge in the terminals 1, 2, ...M as well as in the 
probe. The inequality <^ At, on the other hand, leads 
to that the FCS during At is given by the Levitov-Lesovik 
formula, Eq. Q , with the energy integral replaced by the 
factor dE. To extend the FCS to longer times we write 
all the charges transported during the time interval as 
the difference between the charges in the terminals be- 
fore and after the transfer, Q(t -I- At) — Q{t) = Qi — Qo, 
where Q{t) = {Qi{t), ....,QM{t),Qp{t)). We then make 
use of the composition property of the probabilities 



P (Q2 - Qo) = J dQiP (Q2 - Qi) P (Qi - Qo) 



where we use the notation J dQ = J dQi . . . dQjxjdQp. 
The probability to transfer Q2 — Qo charges during the 
time interval 2At is obtained by multiplying the transfer 
probabilities in each interval At of Q2 — Qi and Qi — 
Qo charges respectively, and integrating over all possible 
values of Qi. 

Extending this to the entire measure- 
ment time r one has P (Qk — Qo) — 

P (Qfe - Qfe-i) n«=i / dQnP (Qn - Qn-i), whcre 
k = T I At. This can be written with the help of Eqs. 
and © as 



p (Qfc - Qo) = y" rfAo n / '^Q" / 



(A2) 



X exp 



■fe-i 



-iA„ • (Q„+i - Q,: 



dEAt 



-ffo(Qri, A„ 



The condition At ^ represents the long time limit, 
when does not depend on time and consequently 
/o^* i?o — AtiJo- The function depends here formally 
on all charges transmitted into the different terminals. 
However, it contains explicitely only the charge on the 
probe via Qp = enp{E)TddE/h, because the distribution 
functions of the terminals 1 to M are externally deter- 
mined, (compare Eq. ([7|l). 

Since the charges Q only change slightly during the 
time interval At Td we can take the continuum limit, 
in which Qn+i — Qn = AtQ„ and the sum in the expo- 
nent in Eq. ljA2|l becomes an integral. We moreover intro- 
duce the standard path integral notation J VQ J VA = 

I rin^i / dQn J rfA„. This gives 



PiQk-Qo) 



X exp [ j dt 



VQVA 



dE 

-iA-Q+—Ho{Qp,A) 



(A3) 



Our interest is the number of charges transported into 
and out of the electronic reservoirs 1 to M, the absorbed 
charge Q° = {Qi, ..,Qm). We can first integrate by parts 
the term in the exponent dtA" • Q" = — /J" dt A" • 
Q° -|- (A^ • Q^ — Ag • Qq) and then functionally integrate 
over Q" to get (5(A"). This functional delta function just 
gives that A" = (Ai, ...Am) are independent of time and 
the functional integrals over A" are reduced to standard 
integrals. Then we have 



P {Qk - Qo) = J dA°exp (-zA'^(Q^ - Qg)) 



(A4) 



X / VQpVXpexp 



dt 



'iXr, 



dE 



HQiQp,A) 



(Al) 



By choosing the boundary conditions so that Ap(t -|- 
t) = we consider the case where the charge Qp is 
not detected, i.e. it is integrated over and we have 
P (Qfc ~Qo) P (Qfe - Qo)- This gives the FCS for 
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the charge in the energy interval dE transported between 
the reservoirs 1 to M during the measurement time t. 

As a last step, we sum up the independent contribu- 
tions from all the energy intervals. Making use of the 
relation Qp — enp{E)TddE/h, this leads to an energy in- 
tegral in the exponent, and a path integral in the field 
np{E,t), and we find the Eqs. (|14|) and H15|) in section 
IV Al describing the generating function 



exp(S') = J 'Dup'DXpexpS with 



(A5) 



S 



dt / dE [-iTdXpfip + Ho] ■ (A6) 



For the voltage probe on instead right from the out- 
set considers the total energy integrated charge Qp — 
(eTd/h) J dEnp(E) and correspondingly the energy inte- 
grated generating function J dEHo. Making use of the 
relation for an equilibrium distribution function np{E) 



we directly arrive at Eq. lEH) 



dVp 



h 



(A7) 



APPENDIX B: SADDLE POINT CORRECTIONS 

In section IV Al the stochastic path integral Eq. p4|l 
was evaluated in saddle point approximation. The func- 
tion S is expanded around the saddle point, see Eq. (fT7|l . 
Taking into account the Gaussian fiuctuations around the 
saddle point, the constraint generating function becomes 
S = 5" + (55, where the upper index stands here for the 
solution at the saddle point. As shown in the following, 
the corrections SS are of order Ts/rd smaller than S'^ and 
can be neglected in the case » considered here. 

For simplicity we consider here energy independent 
scattering and a dephasing probe, other cases are treated 
analogously. In this case, the function S — + 5S can 
be evaluated for each energy interval independently or al- 
ternatively we can directly introduce a factor eV in place 
of the energy integral. We introduce the abbreviations 
A = AAp and n = Arip for the deviations away from the 
saddle point [\^,n^). Following standard path integral 
procedures for saddle point corrections^'^, we insert the 
definition of S*, Eq. H15|) and discretizing time we find for 
the expansion around the saddle point, Eq. (|17|l 



function Hq evaluated at the saddle point. The upper 
bound of the sum in K = r/dt and dt will be taken 
infinitesimal small. 

In order to obtain a purely quadratic matrix 
form it is convenient to use the Fourier trans- 
forms Afc = -^Y.^1-KI2^^^^^ "-^ = 
1 



~7W Y^f=-Ki9 ^ ^'^''^''Ni. One arrives at 



1 eV 

2ir 



J2 (A; Ni A_i N^i) Ti 



l>0 



(B2) 



where the matrix T; has non-zero elements T13 = 
T31 = Adt, T24 = T42 = -Cdt, Ti4 = 

T41 = 2iTd (e^ - 1^ - iBdt, T23 = T32 = 
2iTd (^e^^T^ — 1^ — iBdt. The resulting Gaussian inte- 
gral / dAidNidA-idN-i exp(S' — 5°) is proportional to 
(detT,)"i/2 with 

^/det Ti = {B^ -AC) {dtf +'iTdiTd + Bdt) sin^^ . 

K 

(B3) 

In the limit — > we can expand the sine since only 
small l/K <^ \ contribute in the sum in Eq. (|B2p . 
The term Bdt can be neglected and we arrive at the con- 
straint generating function 



I 



const. 



(B4) 



where k'^ = B^ ~ AC and xi = ^l. The first term 
is the solution 5" at the saddle point which is already 
normalized, 5*° (A = 0) = 0. The second term represents 
the Gaussian corrections. We transform the sum into an 



integral J^i^'^l^'^ + ^fl ~^ 4^ Jo dxln\K'^ + x'^\. This 
contains a divergent contribution which stems from the 
missing normalization. To obtain the correct normalized 
function we take the derivative of SS with respect to the 
counting variables A in the terminals 



dSS 



T dn f°° 2k 
dx- 



T dn f-r, \ 

dA AiTTd dA Jo k'^ + x'^ ~ 4Trf dA 
We integrate again and find for the generating function 



eV 



K 



[-iTd\k{nk+i - nk) 



k=l 





f -A 


iB \ 




nk) 


y iB 


C ) 





(Bl) 



dt 

where the matrix elements A = (^p^), B = (gf^|^) 
and C — ^ ^ contain second order derivatives of the 



S = S° + 5S 



eVr 



HoiXlK) 



4t, 



MA)-k{0)) 



(B6) 

with the correct normalization S{A = 0) = 0. The factor 
r7 ~ Wt^ j^^* ratio between the average time be- 
tween two wave packets = h/eV and the much longer 
delay time of the probe Td, and thus the correction SS is 
small compared to the contribution of the saddle point 
50. 
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APPENDIX C: PHASE AVERAGING AND 
FLUCTUATING PHASES 



APPENDIX D: PROOF OF EQUIVALENCE OF 
ONE-CHANNEL DEPHASING AND PHASE 
AVERAGING 



In section IVIII we compared the result of the dephas- 
ing probe with phase averaging of the generating function 
and found that the two approaches coincide for the case 
of a uniform phase distribution. We pointed out that 
such a phase averaging corresponds to an average over 
ensembles of a chaotic cavities, i.e. an average over static 
disorder. Here we would like to clarify that for the phys- 
ical situation with a phase that fluctuates in time during 
the measurement process, an average of the fluctuating 
phase does in general not correspond to a phase average 
over the cumulant generating function. 

For a phase that fluctuates in time, to find the cor- 
rect averaging procedure it is necessary to compare the 
timescale of the phase fluctuations (say Td) with the mea- 
surement time T. Let us take P{Q) the probability that 
Q charges are transfered during the measurement. In a 
gedanken experiment, one would need to perform m ^ 1 
measurements of the transfered charge to collect enough 
statistics in order to determine the FCS. Three different 
cases can then be distinguished. 

i) Td S> mr. This corresponds to the situation that the 
phase is changing so slowly that it is possible perform m 
measurements without the phase changing. This means 
that all the cumulants, or equivalently the cumulant gen- 
erating function, can be determined for a given phase 
if. Repeating the experiment many times, for different 
phases, one can then obtain the phase averaged cumulant 
generating function {S{ip))^. This is thus equivalent to 
the average over static disorder. 

ii) T ^ Td -€1 mr. In this case the phase is constant 
during a single measurement, however it varies from mea- 
surement to measurement. As a consequence, the charge 
statistics obtained by measuring a large number of times 
is the phase average of the probability distribution itself, 
{P{Q,ip))^, or equivalently the phase averaged moment 
generating function {eyiY>[S {ip)]) ^ . 

iii) Td ^ T. For this situation the phase is fluctuat- 
ing rapidly during a single measurement. The probabil- 
ity distribution giving the statistics outcome of a single 
measurement is independent on the phase, it is already 
given by the phase average {P{Q, (p))ip. Performing a 
large number of measurements thus gives the same result 
as in ii). 

It is important to note that the two averaged functions 
(S)^ and (e^)^ coincide only to linear order in voltage, 
the average of the moment generating function leads to 
modulation contributions in the cumulants which are due 
to the modulation of conductance in time during the ex- 
periment and are of higher order in voltage^^^ 



Here we demonstrate the equivalence between phase 
averaging and the dephasing probe for a single mode 
probe connected to an arbitrary Af-mode scatterer. We 
want to make use of the multi-particle probability expan- 
sion of the cumulant generating function in Eq. (jS)) and 

therefore consider the matrix S^'i^^ . This matrix is formed 
by taking the intersecting matrix elements of the columns 
corresponding to the elements in {a} and the rows cor- 
responding to the elements in {&} from the scattering 
matrix S. For convenience of notation we use Eq. H41|l . 
the scattering matrix U of the entire M + 1 mode system. 
It is then possible to write 

iip 

{0} ^ 

Making use of the general determinant relation 



det 



V W 
X Y 



det{Y) det{V - XY-^W) (D2) 



we can write 
det 



U 



{a} 



u 



Developing this determinant by the last row one obtains 
finally 



det 



det 



«)-e-det(wJ:;,'^/)](D4) 



Taking the modulus square of this determinant and in- 
serting it into the expression for the cumulant generating 
function expanded in terms of multi-particle scattering 
probabilities, Eq. (|7|, we arrive after some algebra at 

H^{<p) = In [A + Bie'^ + B2e-'^] - In [l - Z^/Je'-^lDS) 



Here A is defined in Eq. H26|l above and 

Bi = det [l -I- ni{X\u'^\iU ~ 1) 
B-2. = dot [l + 7ii(A|Zi^AiZ^ - 1) 



(D6) 



with n\ and Ai defined in Eq. (|22|l . We also introduced 
the matrix 



U = 



5oo 
1 



(D7) 
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cal to the result with a dephasing probe, Eq. (|26|l . for a 
single mode probe connected to an arbitrary mesoscopic 
scatterer. 



APPENDIX E: TWO PROBES 



FIG. 10: A four terminal geometry with two probes attached. 
The direction of transport channels is indicated by arrows, the 
intersections marked with yellow bars represent beam split- 
ters. Interference between paths passing through the probes 
in different order as shown on the right is destroyed by de- 
phasing probes, but not by phase averaging. 




FIG. 11: Three barrier setup with two dephasing probes (a) 
or two elastic scatterer (b). While the probes destroy the 
phase memory in each loop individually, the elastic scatterer 
(phase averaging) allow coherent multiple loops. 



Performing the phase average in Eq. H43|) , we arrive at 
the cumulant generating function 



(D8) 



It is then a lengthy but straightforward exercise in deter- 
minant algebra to show that 



B1B2 = B 



(D9) 



where B is defined in Eq. (|26|l . This concludes the proof 
that the result for phase averaging, Eq. HD8|> is identi- 



The agreement between the dephasing probe model 
and phase averaging found for a single mode probe does 
not hold anymore for the case of several probes or of 
probes with more than one transport channel. Differ- 
ences occur already on the level of conductance. A quali- 
tative understanding of the difference between phase av- 
eraging and dephasing probes can be obtained for the 
simplest possible case with two single mode probes. The 
difference arises due to scattering trajectories that visit 
the two probes in different order, but each probe the same 
number of times. These paths give rise to an interference 
term in the conductance that survives phase averaging 
but is destroyed by the dephasing probes. 

The most elementary paths going from say terminal a 
to P, via the probes pi and p2 are then a pi p2 P 
and a p2 Pi (3. An example of a geometry where 
such elementary paths are possible is shown in Fig. 1101 
It contains two dephasing probes pi and p2 or two elastic 
scatterer 0i and (j)2- Transport is unidirectional, and six 
refiection-less beam splitters divide the currents in the 
setup as is indicated by arrows in the figure. The crossing 
in the center of the structure means that the beams are 
not divided but pass one on top of the other. 

The geometry in Fig. 1101 is rather complex. A sim- 
pler three-barrier geometry where phase averaging and 
dephasing give different results were recently proposed 
in Ref. 19 and discussed in section ITXl a schematic pic- 
ture is shown in Fig. ^] Dephasing probes or elastic 
scatterer are connected to the regions between the barri- 
ers. The most elementary paths responsible for the dif- 
ference between phase averaging and dephasing probes 
in this geometry are a — > pi — s- p2 Pi ~^ P2 — > /? and 
a — > pi — > pi — > p2 — > P2 — with a ~ I and /3 = 2 in 
Fig-im These paths are thus more complicated than the 
elementary paths in the geometry in Fig. IIUI We must 
leave it as a further challenge to find additional simple ge- 
ometries for which dephasing and phase averaging leads 
to different results. 
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